Attraction of Spiral Waves by Localized Inhomogeneities with Small- World 

Connections in Excitable Media 
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Trapping and un-trapping of spiral tips in a two-dimensional homogeneous excitable medium 
with local small-world connections is studied by numerical simulation. In a homogeneous medium 
which can be simulated with a lattice of regular neighborhood connections, the spiral wave is in 
the meandering regime. When changing the topology of a small region from regular connections 
to small-world connections, the tip of a spiral waves is attracted by the small-world region, where 
the average path length declines with the introduction of long distant connections. The "trapped" 
phenomenon also occurs in regular lattices where the diffusion coefficient of the small region is 
increased. The above results can be explained by the eikonal equation and the relation between core 
radius and diffusion coefficient. 

PACS numbers: 87.17.-d, 80.40.Ck, 87.17.Aa 



I. INTRODUCTION 



Spiral waves are characteristic structures of excitable 
media that have been observed in many extended sys- 
tems such as reaction-diffusion media [H E 0j aggre- 
gating colonies of slime mold 3, and heart tissues @, 
where they are suspected to play an essential role in car- 
diac arrhythmia and fibrillation. Sudden cardiac death 
resulting from ventricular fibrillation is generated from 
the fragmenting or breakup of spiral waves 0, @. Spi- 
ral waves are prone to a variety of instabilities |QL floL fll| - 
one of which is meander instability jl2l Il3j | , where spi- 
ral tips follow a hypocycloid trajectory instead of moving 
around a small circle. Due to Doppler effect, this spiral 
may undergo a transition from ordered spiral patterns to 
a state of defect-mediated turbulence 11]. In the mean- 
dering regime, the spiral tips can be made to drift and 
controlled by external influences 0| or localized inho- 
mogeneities of defects |l5| . 

After the concept of small-world connections was pro- 
posed by Watts and Strogatz [l(j , it has quickly attracted 
much attention because this kind of connections exists 
commonly in real world, such as in social systems |17). 
neural networks 0] and epidemic problems [l9|. Dif- 
ferent studies show that a little change of the network 
connections can essentially change the features of a given 
medium, and plays a very important role in determining 
the dynamic behavior of a system. 

In numerical simulations, a spatially extended system 
can be approximately regarded as a network consisting 
of a number of sites connected with certain topology. In 
this way localized inhomogeneities can be achieved by 
changing the topology of the network. In heart tissues, 
pacemakers dominate the dynamics of the travelling wave 
behavior and control the heart rhythm. We hypothesize 
that this happens because the characters and the struc- 
tures of the local cells are different from other cardiac 
muscle cells. Could a small-world network describe one 
of the characters of the pacemaker? To answer this ques- 
tion, we changed the widely used regular network in spi- 



ral wave study to a small-world network in part of the 
system to investigate its effects. 

In the next section, we study the effect of a local small- 
world network on the behavior of spiral waves. We show 
that this special region is a dynamic attractor for spiral 
tips. In section III, we compare the effect of the small- 
world network with that of changing diffusion coefficient 
in the local region, and show that they are equivalent. 
We give a discussion and conclude our study in the last 
section. 



II. THE EFFECT OF SMALL- WORLD 
NETWORK 

The model we used is the two-variable FitzHugh- 
Nagumo model [2(| with local nearest-neighbor couplings 
in a region of N\ x N2: 



dui 



dt 



(a - Uij)(u it j - l)l(j,j - v it j + D u V 2 u itj (1) 



dvi 



dt 



= e(buij - v it j) + D V V' 



(2) 



where i — 1, 2, ...,N±, j = 1,2, ...N2; Ujj(t) and Uj,j(t) 
is dimensionless excitable variable and recovery variable 
respectively; D u and D v are diffusion coefficients of the 
two variables. The Laplacien in the last term can be 
approximated as: 

^ 2 Utj = + Ui+i.j + tti,j-i + Uij+i - 4ui,j)(3) 



h 2 
1 



(vi-1,3 + Ui+1,3 + Vi t j-i + Vij+i - 4u ij ).(4) 



When < a < 1, b > 0, D v « D u , e « 1, the 
equation describes an excitable medium, which can be 
regarded as a simplified model for cardiac tissues. In 
the following discussion, we set the control parameters 
as follows: N r = N 2 = 256, a = 0.1, b = 1.0, e = 
0.005, D u — 0.33, D v = 0. No-flux boundary condition 
is applied in the simulation. 
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FIG. 1: Spiral pattern in different media, (a) a regular net- 
work p s = 0; (b) a small-world network, p s — 0.1; (c) spiral 
wave and its tip's motion in the regular network, the white 
curve is the trajectory of the spiral tip; (d) spiral wave and 
its tip's motion in the small-world network. N a = 10. 



Using vertical gradient distribution initial condition, 
we first create spiral waves in regular lattices (Fig. 1(a)). 
In this case, the spiral tip follows a hypocycloid trajec- 
tory, showing a typical sign of meandering state |12j| . 
see Fig. 1(c). We then create a small- world network 
in a small local region O of a size N s ■ N S (N S << iVi), 
where the tip of spiral waves locates. The small-world 
network is created in the following way: With the prob- 
ability p s (0 < p s < 1), we reconnect every edge in the 
region f2 from one of its original vertex to another ver- 
tex chosen randomly in the region [3 (see Fig. 1(b)). 
The change of connections leads to the change of "dif- 
fusion" mode. Supposing node[i][j] is connected with 
node[xi][j/i], node^]^], ■ ■ ■ node [£fc][?/fc], then the "dif- 
fusion" term in the Eqs (1) (2) becomes: 



1 

i 



V 2 Ujj- — -r^(U Xltyi + U X2 ^y 2 + ... + U Xk ^y k k ■ Uij) (5) 



V U M - 12 ( V *l,Vi + V x 2 ,V2+- + Vx k ,y k - k ■ Vij) (6) 

Introducing the small-world network region in the re- 
action medium greatly influences the motion of the spiral 
tip. We find that the small- world network region can at- 
tract the spiral tip as it passes through the region. After 
that, the spiral tip rotates around its boundary, as shown 
in Fig. 1(d). Because the topological structure of the lo- 
cal small- world network is generated randomly under the 
above mentioned rule, the attraction only occurs under 
a certain parameter range with certain probability. To 
characterize the attraction property of the small-world 
network, we use the attraction probability p a as an order 
parameter, which can be obtained by repeating (50 times 
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FIG. 2: (a) transition curves of p a as a function of p s with 
different N s . The line is the sigmoidal fitting p a — 1 — 1/(1 + 
exp((p s — p sc )/dp 3 )); (b) transition curves of p a as a function 
of I with different N 3 . The line is the fitting of p a = 1/(1 + 
exp((/ — l c )/dl)). Here l c is the critical length where p a = 0.5, 
and dl is the transition width. 



in our work) the simulation using the same control pa- 
rameters but with different small- world network connec- 
tions. Our simulation results show that the most influ- 
ential factor to p a is the small- world creation probability 
p s . As show in Fig. 2(a), p a increases with the increase 
of the small- world p s . When p s — (corresponding to a 
regular network), the spiral tip cannot be attracted; when 
p s = 1 (corresponding to a random network), the tip can 
be attracted with probability 1. Between < p s < 1, 
there is a transition where the attraction probability p a 
increases rapidly. The transition point p c can be defined 
as the value of p s when p a = 0.5. 

One of the most important characters of small-world 
network is the reducing of average path length while keep- 
ing the clustering coefficient almost constant 0] . Define 
the normalized average path length I of small-world net- 
work as I — L/Lo, where L is the average path length of 
small- world network in f2 [2l| and Lq is the average path 
length of regular network in Q. I will decrease from 1 to 
e(e > 0) when p s changes from to 1. From numerical 
simulations, we find the same type of transition curve of 
p a as a function of I, as shown in Fig. 2(b), indicating 
that the major effect of the small-world topology to the 
behavior of spiral waves is the decrease of average path 
length. In addition, defining the critical length l c as the 
value of I when p a — 0.5, we find that the critical length 
decreases linearly as the increase of the size of f2, see Fig. 
3(a). 
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III. COMPARE WITH INCREASING THE 
DIFFUSION COEFFICIENT 



Our simulation results suggest that the major effect of 
the small-world network on spiral tip movement comes 
from the long distance connections, which lead to short- 
ening the average path length (/) and increasing the dif- 
fusion speed. If the above suggestion is correct, the phe- 
nomenon of spiral tip attraction should also occur when 
we locally increase the diffusion constant in Eqs. (1) and 
(2). In this part of work, we increase the diffusion coeffi- 
cient in a small circular region Q by D times and keep the 
system with regular connections. In the following discus- 
sion, we will use to denote the diffusion coefficient in 
region fi, and use D® for the region outside of f2, so that 
= D-D^- We find that the spiral tip can be attracted 
and travel around the boundary of region 57 when D is 
large enough and R > Ro.(R is the radius of f2, Rq is 
the core radius of spiral when = D^). At a given 
R, we can define two values D\ and D^: A temporal at- 
traction occurs when D\ < D < £>2- In that case, the 
tip can be trapped for a short period and then escapes; 
the "trapped" time increases with D. When D > £>2, 
the tip can be trapped for a long enough period. Define 
D c as the mean value of D\ and D2 , the plot of D c with 
different R is shown in the Fig. 4. 

The above simulation results suggest that the increase 
of diffusion speed in the small region is responsible for the 
attraction of the spiral tip. To quantitatively compare 
the two systems, we analyze the diffusion terms of the 
two systems. In the small-world network, because of the 
long distance connections, the average distance between 
nodes declines as p s increases from to 1. In a network 
model of a reaction-diffusion system, this effect can be 
in a sense translated from the decrease of the average 
path length between nodes while keeping the distance of 
two neighboring nodes h constant, to the decrease of the 
step length h while keeping the network regular. The 
normalized average path length / can also describe the 
relative change of h. From this argument, the diffusion 
items of the small-world network can be expressed as: 



D 



:(Vi-l,j + Vi+lj + + V itj+1 - AVij) (8) 



v (lh)- 



At the critical point, for a fixed R, the diffusion terms 
in two systems should be the same. So that ■ 1/h 2 — 
D c -D°-l/h 2 = D° u -l/(l c -h) 2 , which gives D c -l 2 c = 1 . As 
presented in Fig. 3(b), our simulation fits this analysis 
within the range of error. This result indicates that our 
proposition is reasonable. The attracting effect of the 
small-world network comes from the decline of I inside 
the inhomogeneous area f2. 




FIG. 3: (a) The critical length l c as a function of the region's 
radius R, (R = Ns/2), the error bar is dl,(see the caption of 
Fig. 2(b) about the fitting); (b) D c ■ 1% as a function of R. 



IV. DISCUSSION 

A question should be answered before fully under- 
standing the effect of the small-world network in the dy- 
namics of spiral tips: what is the mechanism for the spi- 
ral tip attraction? In the following discussion, we give an 
explanation with the eikonal equation, Luther equation 
|22j , and the relation between diffusion coefficient and of 
spiral core radius. According to the analysis of the spi- 
ral tip dynamics given by Hakim and Karma |2^| , for the 
steady rotational movement of a spiral tip in an excitable 
medium, the core radius R as a function of diffusion co- 
efficient can be written as: 



R=^(- 
co B, 



b-K 



ID, 
W 



3/2 



(9) 



Where D u is the diffusion coefficient of activator (it); cq 
is the speed of plane wave; b, K and B c are all constants. 
W is the constant width of the excited region. In the 
simulation, We assume that at the boundary of Q exists 
a "virtual" gradient region of D u which links the outside 
and inside regions. For a given R of the region fl(R > 
Rq = R(D^)), the "trapped" motion of the spiral tip 
requires a specific value of D Ul satisfying equation (jSJ. 
When D u < D^, the spiral tip will enter the gradient 
region where the system can find the required D Ul so that 
the spiral tip will rotate around the gradient area at the 
boundary of SI; on the other hand, when D u > D^, the 
"trapped" motion can not be sustained by the central 
region. From this argument, at critical point, we will 
have D U (R) = D% = D c ■ D°. As shown in Fig. 4, Our 
simulation results are consistent with this analysis with 
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In our system, assuming there exists a continuous 
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FIG. 4: The critical diffusion coefficient D c as a function of 
R, where the line is the best fitting with Equation The 
error bars are estimated using D\ and L>2. 



the range of error. 

To prove the "trapped" state of tip motion is a sta- 
ble state, we apply the eikonal equation, which deter- 
mines the relation between the curvature of a travel- 
ling wave front and its speed in an excitable medium, 
and the Luther relation, which describes the relation 
between the speed of chemical waves and the diffusion 
coefficient of activator j2^. the eikonal equation is: 
N = C — D u ■ k, where N is the normal wave speed, 
K is the local curvature of the wave front; the Luther 
equation is: c = a ■ ^/D u , where a is a constant. Insert 
the Luther equation into eikonal relation we have: 

N = a- y/D^-Dv-K (10) 
Taking partial derivative of R in Eq. (|10f> . we get: 



ON 
~dR 



A 



a (A, 



-1/2 



dD u 



l 2 -v QR 

At the spiral tip we have ./V = 0, so that: 

1 dD n 



ON 
9R \tip 



' 2 K tip 



dR 



(11) 



(12) 



dDu 
dR 



< 



change of D u at the boundary of region ft, we have 
0. Thus the equation ijT^} indicates that \ tip ) 0. 
That means, if we introduce a small deviation from the 
"trapped" motion of the spiral tip, the system will re- 
turn to the "trapped" state spontaneously, because we 
have N < inside the region O, and N > outside the 
region Q. (fi points to the center of Q region). 

In conclusion, we find that in an excitable system lo- 
cal change of topological structure can trap the spiral tip. 
This ability comes from the increase of diffusion speed. 
We prove this by increasing the diffusion coefficient. We 
also give a theoretical explanation using the eikonal equa- 
tion, Luther equation, and the relation between core ra- 
dius and diffusion coefficient, which fits well with the re- 
sults of simulation. We should note that there are other 
situations where the tip of spiral waves can be trapped 
in a given area. For example, Lazar et al. reported that 
self-sustained chemical waves can rotate around a cen- 
tral obstacle in an annular 2-D excitable system, and the 
wave fronts in the case of an annular excitable region 
are purely involutes of the central obstacle in the asymp- 
totic state 24] . Obviously this phenomenon is beyond 
our analysis. So that more work should be done to fully 
understand the attract effect of local inhomogeneities in 
an excitable reaction-diffusion system. 
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